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Methods and App aratus of Sample Analysis 

The present invention relates to methods and apparatus for the analysis 
of a sample. It relates particularly to x-ray analysis, and has applications in both 
' 5 the medical and non-medical fields. It may be used in the analysis of organic 
and inorganic substances, and crystalline and amorphous materials. 

X-ray analysis in the form of x-ray radiography is well-known in the 
medical field, and is a widely used technique for medical imaging. 

In its simplest form, a radiographic device comprises an x-ray source and 

10 an x-ray sensitive film arranged opposite to one another on either side of a 
patient. X-rays from the source are directed at the patient's body; and those 
that pass through are recorded on the film. Since bone and soft tissue absorb 
different amounts of radiation, an image of lighter and darker areas is captured 
on the film, and can be used to identify internal structure of the patient. 

1 5 A problem with this technique is that organic materials are generally 

weak absorbers of X-ray radiation, and so distinguishing between different 
types of soft tissue can be difficult. For example, it can be difficult to distinguish 
between healthy soft tissue and a tumour, especially in a tumour's early stages. 
X-ray techniques based on phase imaging have recently been promoted 

20 for overcoming this problem. Rather than use x-ray absorption, these methods 
attempt to provide better image contrast by utilising instead the phase shift of x- 
rays as they pass through regions of the sample of different refractive index. 
Phase imaging techniques can provide improvement over radiographic 
techniques, and can allow for better definition of soft tissue. 

25 A problem with both radiographic and phase imaging techniques is that 

in order to use relatively short exposure times, they generally require the 
incident radiation to be of a relatively high intensity. Damage of the object 
under examination by the x-rays is therefore of concern. This is especially so in 
medical imaging, where strict protocols must be observed as to radiation 

30 dosages and the like. 

A further point of note is that it can be difficult and expensive to obtain 
high resolution in both techniques. For example, even the use of the most 
expensive CCD cameras (CCD refers to charge couple device) on the market 
today still imposes undesirable limits on image resolution and dynamic range. 
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The present invention aims to provide an alternative x-ray analysis 
technique that, in its various aspects and embodiments, can provide a number 
of advantages both in the medical and non-medical fields. 

Viewed from a first aspect, the present invention provides a method of 
5 analysis of an object, the method including the steps of: 

generating non-planar penetrating radiation; 

diffracting the radiation from a monochromator to provide a beam of 
monochromatic penetrating radiation; 

irradiating a portion of the object with the beam; 
1 0 diffracting radiation that passes through the object onto a detector from 

an analyser; 

rotating the analyser through a plurality of angular positions; and 
measuring the intensity of the radiation incident on the detector as a 

function of analyser position. 
1 5 Viewed from a second aspect, the present invention provides apparatus 

for the analysis of an object, the apparatus including: 
a source of non-planar penetrating radiation; 

a monochromator for diffracting the non-planar penetrating radiation to 
provide a beam of monochromatic penetrating radiation; 
20 a detector for detecting radiation that passes through the object; 

an analyser for diffracting radiation that passes through the object onto 
the detector; 

means for rotating the analyser between a plurality of angular positions; 

and 

25 means for recording the intensity of the radiation incident on the detector 

as a function of analyser position. 

Each intensity value detected through the present invention (i.e. for each 
angular position of the analyser) corresponds to the intensity of the x-ray 
photons that are diffracted by the analyser onto the detector from across the 

30 whole of the radiation beam that has penetrated the sample object. The 
resultant intensity profile includes information as to both the absorption and 
refraction of x-rays by the sample, and facilitates the determination of for 
example the shape and structure of the analysed portion in a single experiment. 
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As is preferred, the present invention allows for the determination of a 
complex scattering function of the analysed portion from the intensity 
measurements. This can provide the complex refractive index of the portion of 
the object under analysis as a function of cross-section transverse to the 
5 direction of propagation of the beam. 

The invention contrasts with the prior art radiographic and phase contrast 
techniques that generally view the object under analysis either in a purely 
absorptive or in a purely phase manner. 

The present invention is further distinct from the prior art in that it can be 
10 thought of as obtaining and analysing information in "reciprocal space" rather 
than in "real space". 

The present invention facilitates the analysis of an object at high 
resolutions in an inexpensive manner. For example, there is no need to use an 
expensive high resolution CCD detector, as the detector used in the present 
15 technique need not provide spatial resolution (spatial resolution is instead 
provided by the rotation of the analyser, and may be determined by the full 
aperture of angular rotation). 

The present technique can provide higher resolution than the prior art, 
and may for example provide resolutions in the region of 0.05 - 0.5 nm. 
20 The present invention is also able to efficiently utilise the x-ray photons 

emerging from the portion of the object analysed, since it effectively integrates 
the emerging x-rays across the whole width of the analysing beam, and the 
detector may start counting from a single photon of light. The present invention 
thus facilitates the use of x-ray radiation of lower intensity, and may provide 
25 better data statistics and so analysis quality than imaging techniques. 

It should be noted that the present invention analyses non-Bragg 
diffracted x-rays exiting from the sample object, and thus is distinct from e.g. 
standard crystallography and the like. 

Generally, an analysis of non-Bragg diffracted waves would be 
30 considered of no use, as it was not thought possible to obtain any meaningful 
information from such an analysis. Indeed, the idea of obtaining diffraction 
effects at all for non-periodic structures, such as amorphous materials was not 
entertained. 
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The inventor of the present analysis apparatus and method has not only 
found that meaningful information can be obtained from a study of non-Bragg 
diffracted waves from both crystalline and non-crystalline samples, but has also 
found that quantitative information is achievable. 
5 Further, not only has the inventor found this, but also has determined 

contrary to what might be thought, that such information is obtainable even 
when using a standard laboratory source (as opposed for example to when 
using beams of ultra-high intensity and ultra-low divergent radiation that are 
only found in the most expensive of synchrotrons). 
10 The inventor has thus been able to provide an analysis technique that 

offers a practicable and high-resolution alternative to the prior art imaging 
methods. 

The use of laboratory sources (as opposed to top-end e.g. 1 kilometre 
long synchrotron radiation sources) may be characterised by the generation of 

15 non-planar penetrating radiation, which is then collimated by the 

monochromator (with a resultant loss in beam intensity). Non-planar 
penetrating radiation in the present application may be considered to be 
radiation having an angular divergence of greater than 10 seconds of arc when 
incident at the monochromator. 

20 In a preferred embodiment, the x-ray source is a characteristic line 

source such as an x-ray tube. In such sources, material is bombarded by 
electrons to generate the x-rays. The material may be e.g. copper, 
molybdenum, silver or tungsten. 

In an especially preferred embodiment, the source is a rotating anode 

25 source. Such sources are similar to standard x-ray tube sources, but the 
material bombarded is rotated. Rotating anode sources are able to provide 
relatively high intensity beams so as to provide greater resolution and better 
signal to noise ratios in the intensity data. 

For comparison, an x-ray tube typically generates x-rays having an 

30 energy range of from about 5 to 1 50 keV and having an intensity of about 1 0 6 
photons/sec, whilst a rotating anode source typically generates x-rays having an 
energy range of from about 5 to 60 keV and having an intensity of about 10 8 
photons/sec. Synchrotron radiation typically has an energy range of from about 
1 to 100 keV and an intensity of about 10 12 " 18 photons/sec. 
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Other laboratory sources of x-rays could also be used. These could 
include for example small-scale synchrotrons, e.g. which are able to fit into a 
room of a building, but which, because of their small beam lengths, may still 
require optical elements, such as a monochromator, to provide a beam of low 
5 angular divergence. 

The analyser may take any suitable form, and for example may comprise 
a crystal oriented so that only x-rays of a set angle of incidence or range of 
angles and desired wavelength are diffracted towards the detector. It may 
comprise a single or double crystal that may be made of e.g. silicon or 
10 germanium. Silicon is preferred, and provides a stable and effective analyser. 
The crystal analyser is the analyser of choice, but alternatives are 
possible. For example, the analyser could be in the form of an x-ray mirror. 

The monochromator may also take any suitable form, and may take a 
form similar to that of the crystal analyser, e.g. a single or double crystal of 
1 5 silicon or germanium, or an x-ray mirror. 

The monochromator and analyser pair provides the system with a set 
optical resolution, determined from e.g. the angular divergence of the beam 
leaving the monochromator and e.g. the angular acceptance of the analyser. 
The optical resolution will generally be the angular divergence of the beam after 
20 passing through the monochromator, or the angular acceptance of the analyzer, 
whichever is greater. 

It is preferred that the monochromator and the analyser are of the same 
type and have equivalent angular divergence and angular acceptance values. 
However this is not essential for the operation of the invention. 
25 The actual optical resolution of the apparatus can be taken from a 

rocking curve of the apparatus, i.e. an intensity profile of detected intensity 
versus analyser position for the apparatus without a sample in place. The 
optical resolution may be defined by the Full Width at Half Maximum of the 
rocking curve (FWHM). 
30 The monochromator and analyser pair may be arranged so that the 

analyser is at its zero angle position when it is diffracting x-rays into the detector 
at the Bragg angle. 

Preferably, x-rays pass through a slit prior to their incidence on the 
object, the slit size A (in a direction transverse to the direction of propagation of 
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the beam and in the direction in which the refractive index is being measured) 
being such that: 

A<A/89, 

where A, is the wavelength of the incident radiation, and 59 is the 
5 optical resolution of the apparatus used (e.g. the monochromator/analyser pair 
resolution, e.g. determined from the FWHM of the rocking curve). 

A point to note here is that the optical resolution 59 is determined by the 
apparatus, rather than being a quality of the beam itself (which would be the 
case if a top-end synchrotron beam where used). 
10 This use of a slit width within the above constraint ensures that the 

intensity profile detected is able to provide analytical functions for determining 
the object's refractive index, and so allows a quantitative analysis of the cross- 
sectional refractive index. 

Typical slits widths may range from e.g. 50 to about 500 microns. A 
15 copper x-ray tube may for example utilise slit widths of between about 50 and 
about 300 microns, whilst a molybdenum source may utilise slit widths between 
about 50 and about 500 microns, and a silver source slit widths of between 
about 50 and about 300 microns. Preferably, the slit width is greater than 50 
microns, and more preferably greater than 100 microns. It is preferably less 
20 than 500 microns, and may for example fall within the range of about 250 to 
about 300 microns. 

The slit member that defines the slit is preferably made from tantalum, 
and the edges of the slit are preferably flat and parallel to the direction of 
propagation of the beam. This avoids the possibility that x-rays will penetrate 
25 the slit member about the edges of the slit, and cause undesired diffraction 
effects in the measured intensity profile. 

As well as providing quantitative information, the intensity profile may 
merely be plotted e.g. as Log(intensity) against angular position of the analyser, 
and the resulting plot may be reviewed qualitatively e.g. against a similar plot for 
30 a reference sample. 

The size of the slit in the second transverse direction may take any 
suitable value, and may be e.g. between about 0.1 mm and about 10 mm. The 
particular size may be dependent on the sample under analysis, and any 



WO 03/106983 PCT/AU03/00748 

7 

desired spatial resolution in the second transverse direction (e.g. if an intensity 
profile were to be obtained in this direction also). 

Preferably the crystal analyser is rotated in a plurality of incremental 
steps, each step being through an angle a : 
5 a < 69/2 

where 89 is again the optical resolution of the apparatus. 
This constraint allows the complete determination of the complex 
scattering function. It takes into account the need to measure 2N points of a 
real function (i.e. intensity) in order to obtain N points of a complex function 
10 (which has both real and imaginary components). 

A typical step may be e.g. 0.05 or 0.01 arc.secs, although other step 
sizes are equally possible. 

The detector may take any suitable form, and, as said, need not provide 
any spatial resolution. Preferably, the detector comprises a PIN (P-intrinsic-N) 
15 diode detector. PIN diode detectors are able to provide a linear response over 
a large dynamic range. This can be an advantage in the present technique, as 
the scattered intensity profile to be recorded typically has a large variation in 
intensity, with useful information being contained in both high and low intensity 
areas: A typical profile of the logarithm of intensity against analyser angle 
20 comprises a large central peak with small but significant variations in intensity in 
the tails of the peak. As an alternative, a scintillation counter or the like may be 
used. 

Preferably, a complex scattering amplitude of the irradiated portion of the 
object (in reciprocal space) is calculated from the detected intensities, and a 
25 complex scattering function of the irradiated portion is determined by taking an 
inverse Fourier Transform of the complex scattering amplitude. 

The method preferably includes the steps of: normalising the detected 
intensities; calculating the modulus of the complex scattering amplitude from the 
normalised intensity; calculating phase information (e.g. a minimal phase) of the 
30 complex scattering amplitude from the modulus of the complex scattering 

amplitude; and determining the complex scattering amplitude from the modulus 
and the phase information. 

The analysis of the detected intensity profiles may be carried out by any 
suitable means, and may be implemented in software running on any suitable 



WO 03/106983 PCT/AU03/00748 

8 

computing apparatus, for example a personal computer, as would be well 
understood by a person skilled in the art. The computing apparatus could 
control the whole analysis operation and co-ordinate control of the analyser and 
detector, e.g. control the rotation of the analyser and the recordal of the 
5 detected intensities against angular position of the analyser. 

The sample may be positioned relative to the slit, such that significant 
variation in the sample refractive index is expected only in the diffraction plane 
of the monochromator-analyser pair. 

In one preferred embodiment, the object is translated relative to the 
10 beam, e.g. in a linear fashion, between obtaining profiles, so as to provide an 
analysis of the object over a number of beam widths. The results of the 
separate analyses of the various portions of the object may then be combined 
so as to provide a complex refractive index profile of the object over a large 
cross-sectional area. An object can therefore be mapped over a large area at a 
15 high resolution. 

The present invention provides a refractive index profile for the sampled 
portion of the object in a direction transverse to the direction of beam 
propagation, the resulting profile also being the refractive index integrated over 
the distance that the beam propagates through the sample. In order to provide 
20 three-dimensional information, tomographic techniques may be employed. 

It should be noted that the restriction of the slit width in accordance with 
the constraint discussed above is in itself an important feature of the present 
invention, and provides an intensity profile that is analytical and so allows for a 
quantitative analysis of the refractive index of the sample. 
25 Accordingly, viewed from a further aspect, the present invention provides 

a method of analysis of an object, the method including the steps of: 
generating penetrating radiation; 

diffracting the radiation from a monochromator to provide a beam of 
monochromatic penetrating radiation; 
30 passing the beam of radiation through a slit, the slit size A (in a direction 

transverse to the direction of propagation of the beam) being calculated such 
that: 

A<M59 
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where X is the wavelength of the incident radiation, and 80 is the optical 
resolution of the apparatus used in the method; 

irradiating a portion of the object with the beam; 

diffracting radiation that passes through the object onto a detector from 
5 an analyser; 

rotating the analyser through a plurality of angular positions; and 

measuring the intensity of the radiation incident on the detector as a 
function of analyser position. 

The present invention may also be seen as providing an apparatus for 
10 the analysis of an object, the apparatus including: 

a source of penetrating radiation; 

a monochromator for diffracting the penetrating radiation to provide a 
beam of monochromatic penetrating radiation; 

a slit member defining a slit through which the beam passes prior to the 
15 beam's incidence on the object, the slit size A (in a direction transverse to the 
direction of propagation of the beam) being calculated such that: 

Azme 

where A, is the wavelength of the incident radiation, and 80 is the optical 
resolution of the apparatus; 
20 a detector for detecting radiation that passes through the object; 

an analyser for diffracting radiation that passes through the object onto 
the detector; 

means for rotating the analyser between a plurality of angular positions; 

and 

25 means for recording the intensity of the radiation incident on the detector 

as a function of analyser position. 

Viewed from a still further aspect, the present invention provides a 
method of analysis of an object, the method including the steps of: 

irradiating a portion of the object with a beam of monochromatic x-rays; 
30 detecting the intensity profile of an angular spectrum of the x-rays 

emerging from the irradiated portion; and 

determining a complex scattering function for the irradiated portion of the 
object under analysis. 
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Viewed from another aspect, the present invention provides a method of 
analysis of an object, the method including the steps of: 

irradiating a portion of the object with a beam of monochromatic x-ray 
radiation; 

5 diffracting x-rays emerging from the sample into an x-ray detector using 

an analyser means; and 

obtaining an angular spectrum of non-Bragg diffracted x-ray intensities 
as a function of angular position of the analyser means. 

Viewed from a further aspect, the present invention provides a method of 
10 analysis of an object, the method including the step of collecting generic (non- 
Bragg diffracted) x-ray diffraction data from a portion of the object and analysing 
the data to obtain a complex refractive index of the sampled portion in a 
direction transverse to the beam propagation. 

These further aspects may extend to apparatus for applying the 
15 methods, and may include any of the features of the preceding aspects of the 
present invention. 

Embodiments of the present invention will now be described, by way of 
example only, with reference to the accompanying drawings. It is to be 
understood that the particularity of the drawings does not supersede the 
20 generality of the preceding description of the present invention. 
In the drawings: 

Figure 1 is a schematic diagram of x-ray analysis apparatus in 
accordance with a first embodiment of the present invention; 

Figure 2 shows scattered intensity profiles for three analysed samples as 
25 a graph of Log(intensity) against detection angle; 

Figure 3 shows the scattered intensity profile for a sample when placed 
at two different positions relative to a beam slit; and 

Figure 4 represents the profile of the real component of the complex 
refractive index for the two profiles of Fig. 3. 
30 Referring to Fig. 1, an object 1 to be analysed, e.g. a blood vessel, will 

have regions of differing refractive index, e.g. an outer region 2 consisting of the 
blood vessel wall, and an inner region 3 consisting of the blood itself. 

Overall, in order to analyse the internal structure of the object 1 , a beam 
of x-ray radiation (generally labelled 4) is passed through the object 1, and the 
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emerging radiation that is within the acceptance angle of a crystal analyser 5 is 

diffracted onto an x-ray detector 6. The detector 6 records an intensity profile of 

this radiation as a function of angular position of the crystal analyser 5, and 

suitable control means 7, such as a computer running suitable software, 
5 analyses the resulting profile to provide a complex refractive index profile for the 

object 1 across the width of the beam 4 in the x-axis direction (transverse to the 

direction of propagation of the beam 4). 

In more detail, in the apparatus shown, an x-ray source 8, such as an Mo 

x-ray tube with Kq = 0.709nm, generates an x-ray beam 4a that is diffracted 
10 from a crystal monochromator 9 and passed through a slit member 10 so as to 

provide a narrow collimated beam of low divergence monochromatic x-rays 4b 

that is incident on the sample object 1 . 

The monochromator 6 is an Si crystal providing a (333) asymmetric 

reflection, and the crystal is cut such that (111) atomic planes are at 1 8.5° to the 
15 surface. 

The slit member 10 is made of e.g. tantalum and the sides of the slit are 
flat and parallel to the direction of propagation of the beam so as to prevent x- 
rays from passing though the edges of the slit member 10 about the slit opening 
and causing undesirable diffraction effects. 
20 The slit size, A, is selected to satisfy the relationship: 

A < A/56 

where Ms the wavelength of the incident radiation, and 86 is the 
optical resolution of the apparatus as determined by e.g. the angular divergence 
of the beam after the monochromator or the angular acceptance of the 
25 analyzer. In the present case, the optical resolution may be determined by the 
Full Width at Half Maximum (FWHM) of the rocking curve of the slit member 10, 
the rocking curve being the intensity profile obtained by the apparatus in the 
absence of a sample to be analysed. 

The x-ray beam 4b is absorbed and scattered in its passage through the 
30 sample 1 , and the emerging x-rays 4c that fall within the acceptance angle of 
the crystal analyser 5 are diffracted towards the detector 6. 

The crystal analyser 5 may be of the same configuration as the 
monochromator 9, and may comprise an Si crystal providing a (333) 
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asymmetric refraction, the crystal being cut such that (111) atomic planes are at 
18.5° to the surface. 

The crystal analyser 5 is arranged opposite to the entry point of the x- 
rays into the object 1, and is mounted to rotate in incremental steps about an 
5 axis 5a over a set range of angular positions so as to present itself at a plurality 
of different angles to the x-ray beam 4. 

The means 1 1 for rotating the crystal analyser 5 may comprise e.g. a 
goniometer, and may be controlled by the control means 7, which also records 
the intensity of x-rays received at the detector 6 as a function of the angle of the 

10 analyser 5 from its zero position. The "zero" position of the analyser 5 relative 
to the beam 4 corresponds to the Bragg reflection angle of the analyser crystal, 
and in this position, the analyser is in this embodiment substantially parallel to 
the monochromator 9. 

In use, the crystal analyser 5 is rotated between e.g. -30 arc.secs and 

15 +30 arc.secs of the Bragg reflection angle, in steps of e.g. 0.05 arc.secs. The 
steps should be less than half the optical resolution of the apparatus, e.g. as 
defined by the monochromator and analyser pair (This takes into account the 
need to acquire 2N real values in order to compute N complex values). 

At each step, an intensity reading is recorded. This may be for a set 

20 period of time, for instance between about 1 to about 1 0 seconds. Alternatively, 
the intensity of the radiation at each step may be measured until a set number 
of photons is reached, e.g. 1000 photons. In such a case it may take e.g. from 
about 10 to 20 seconds to measure an intensity value at an extreme angular 
position, and e.g. from about 0.1 to 0.5 seconds to measure intensity at a point 

25 close to the peak of the profile. Each resulting intensity value may be modified 
to take account of the time taken to obtain it, e.g. by dividing the intensity value 
by the measurement time. 

The detector 6 may be a PIN diode detector, which is able to provide a 
linear response over a large dynamic range. This is useful for the present 

30 technique, as the detector 6 needs to record both a large central peak of 

substantially unabsorbed and unrefracted radiation, as well as meaningful small 
perturbations in the tails of the profile caused by phase shifts in the x-rays as 
they pass through regions of the object of different refractive index. Alternative 
detectors may also be used, e.g. a scintillation counter. 
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Examples of intensity profiles that may be obtained through the use of 
the apparatus of Fig. 1 are shown in Fig.2, which plots Log (Intensity) against 
crystal analyser angle. 

The profiles A-C are of a polyethylene tube (approximately 6 mm in 
5 diameter and sealed at both ends) filled with air, water and butter respectively, 
so as to roughly simulate a blood vessel configuration: The air is used to 
imitate an empty blood vessel, whilst the water and butter are used to imitate a 
blood vessel filled with clear blood and organic fat respectively. 

The profile RC provides an analyser rocking curve corresponding to no 
1 0 sample. The rocking curve may be used to calculate the optical resolution of 
the apparatus by determining the FWHM. 

It should be noted that the scans are offset by about 1 decade for clarity, 
and that the sample may be immersed in water to reduce the change in 
refractive index that would otherwise occur at the interface of the sample object 
15 with air. 

Inspection of such plots can provide qualitative information as to the 
object under analysis, e.g. one may be able to tell from such a plot if a blood 
vessel is abnormal. 

The present invention can however also provide quantitative information 
20 as to the structure of the object from a suitable analysis of the intensity profile. 

Overall, the concept is to obtain (in reciprocal space) the complex 
diffraction amplitude of the portion of the object analysed by the x-ray beam 4, 
and from this to obtain the complex refractive index profile of the analysed 
portion across the width of the beam. This complex refractive index will have 
25 information on both the absorptive and refractive features of the object under 
analysis, and so will provide information on the structure of the analysed portion 
of the object. 

The analysis relies on the analytical properties of the complex diffraction 
amplitude, which is provided by the use of a slit width A within the above-noted 
30 constraints. 

The refractive index profile obtained will be that along the x-axis of the 
portion of the object that is analysed (integrated over the length of the object 
through which the beam passes),, and assumes that the object is homogeneous 
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or has a slowly varying refractive index in the y-axis direction (In Fig. 1, this 
direction is out of the plane of the paper). 

As a specific example, a profile may be analysed by firstly normalising 
the measured intensity, l(8), by the square modulus of the complex scattering 
5 vector: 

Q 2 =q 2 +M 2 , 

where q= 8/A and p is the linear attenuation coefficient of the 

material. 

The normalised intensity is: 

10 T(Q)^l(Q)x0 2 . 

Next, the modulus \R(Q)\ of the complex diffraction amplitude, R(Q) , is 

calculated as the square root of the normalized intensity: 



R(Q) = yl7(Q)- 

The minimal phase is then calculated from the square root of the 
1 5 measured intensity using a logarithmic Hilbert transform. This can be 

implemented using the relation between Hilbert and Fourier transforms (see e.g. 
R.N. Bracewell, The Fourier transform and its applications (McGraw-Hill, New 
York, 1986): 

^(Q) = -lpjH^, 

20 where P is the Cauchy principle value of the integral. 

The solution for the Complex Diffraction Amplitude can then be 
synthesised as: 

Next, the inverse Fourier transform of R(Q) is taken: 

The complex refractive index profile of the sample across the width of the 
beam can then be obtained by taking the integral: 

A 

X(*)=S%(x)dx. 

0 

An unambiguous solution to the complex refraction index profile may be 
30 obtained either by using intensity profiles collected for two different x-ray 
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energies or by using a single intensity profile compared to previously obtained 
benchmark profiles for similar samples. The results obtained for the two profiles 
can then be compared so as to identify physically real roots from artefacts of the 
computation. 

5 In one experiment, the apparatus of Fig. 1 was used to analyse a sample 

of nylon (C 6 Hi3N0 2 , density - 0.4336 g/cm 3 ) of 80 micron diameter. 

The sample was arranged to occlude the slit in the slit member 10 such 
that the whole sample is "seen" through the slit plus a bit of air. 

The apparatus had an optical resolution of 80 =0.12 arc.sec. = 5.87x1 0" 7 
10 rad, and the slit aperture was selected to be 1 10 microns, which satisfies the 
condition that the slit aperture must be smaller than: 

A = )Jd9 = 7.09 xlCT 11 / 5.87 xKT 7 =120 micron. 
The angular step in the diffraction pattern, i.e. the angular step of the 
crystal analyser 5 was chosen to be 0.05 arc.sec. - at least two times smaller 
15 than the optical resolution. 

The intensity profile (generic optical diffraction pattern), 1(9), was 
determined from the sample by the angular scanning of the analyser from -10 
arc.sec. to +10 arc.sec (giving 401 experimental points). This gave a total 
angular scan aperture: 
20 Qapenure = 20 arc. sec. = 9.7 x 1 0" 5 rad. 

This corresponds to a resolution for the reconstructed profile of the 
refractive index of the sample across the beam width of: 

Ax = A/6 = 7 x 1 0" 11 / 9.7 x 1 0~ 5 = 0.72 micron. 
Further examples of intensity profiles obtained using the apparatus of 
25 Fig. 1 are shown in Fig. 3. They are for a model blood vessel comprised of a 1 
mm diameter grass straw cylinder located partially in the x-ray beam 4 (only one 
edge of the straw was within the beam, the other edge was blocked by the slit 
member 10). 

The sample was immersed in water, so that the total thickness of the 
30 sample cell was about 1 .5 mm, with 2 x 50 micron Kapton™ windows to the 
sample cell. 

The x-ray radiation was Mo Kqi = 17.5 keV. 
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The monochromator/analyser pair used a single silicon dislocation-free 
crystal with an asymmetric cut that had an asymmetry factor b = 0.038. Bragg 
reflection was Si 333, and Darwin's width was 0.12 arc.sec. 

The slit size was 0.1 mm (x-axis, A dimension) x 5 mm (y-axis). 
5 The two profiles A' and B' shown in Fig. 3 are for two positions of the 

sample relative to the 100 wide slit. Profile RC is a rocking curve for the 
apparatus corresponding to no sample, i.e. an open slit. 

In scan A\ the sample occludes the slit by 58 ± 2 pm, and in scan B\ the 
sample occludes the slit by 40 ± 2 pm. 
1 0 Fig. 4 shows the reconstructed profile of the real component of the 

complex refractive index of the sample (proportional to the thickness/density 
profile in the direction of the incident x-rays). 

The drop on the right side of the profile corresponds to the slit edge, and 
the drop on the left-side corresponds to the edge of the sample. The left-side 
15 drop has shifted by about 15 micron between the two profiles, corresponding to 
the movement of the sample between measurements of the profiles. 

The spatial resolution of the reconstruction was 1 .5 pm. 

The basis for the present approach to x-ray analysis is as follows. In the 
case of a non-crystalline sample, the complex scattering function of the sample 
20 may be expressed as: 



where n is the map of the complex refractive indices, n(x,y,z) = 1- 
d(x,y,z)-ifi(x,y,z), and z is the direction of the incident wave propagation, for an 
ideally monochromatic source of wavelength A. 



detector 6) is the square of the modulus of the complex scattering amplitude, / = 
IT(Q)! 2 , where T(Q) is given as the Fourier transform of the complex scattering 
function t(x): 




25 



The experimentally measured scattered intensity (i.e. that obtained at 



T(Q) = \T(Q)\ expOWQ)) = j t(x) exp(2mQx) dx 



all space 



30 



where T(Q)=u(q r ,qi)+iv(q r ,qi) and q r ,qi are the real and imaginary 
parts respectively of the complex scattering vector Q=qr+iq. 
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The inverse problem of determining the complex scattering function t(x) 
from the experimentally observed scattered intensity profile l(Q) depends on the 
theoretical approach to the relation between the modulation function and the 
observed intensity function. 
5 The phase q>(Q) of the experimentally observed x-ray scattering profile 

may be retrieved via a logarithmic dispersion relation, an approach that is valid 
under the kinematical theory of x-ray scattering: 

where OF (m=0,1 ,2,... M-1) are the zeros of IT(Q)/ in the upper 
10 half of the complex plane, and P is the Cauchy principal value of the integral. 
The zeros of the scattered amplitude, QT, are of unknown number, M, and may 
in principle be infinite in number. However, the analyticity of the CD A allows its 
unique representation as a complex polynomial function: 

15 in which the number of zeros is limited to the discrete number of points K 

in the experimentally collected data set. 

Once the x-ray phase, <p(Q) 3 is determined, e.g. the minimal phase 
change term, an inverse Fourier transform of the calculated complex scattering 
amplitude T(Q) can be used to find the complex scattering function t(x) and 

20 provide information on the thickness/composition profile. 

Overall, the present invention can provide a method and apparatus for 
the x-ray analysis of a sample, which is able to provide high resolution and good 
contrast without great expense, and can provide information as to both 
absorptive and refractive features of the sample in a single experiment. The 

25 invention utilises an angular spectrum of non-Bragg diffracted x-rays obtained 
from standard laboratory x-ray sources, and is not limited for example to 
expensive top-end synchrotron sources or the like. The invention can ensure 
analyticity of the resulting profiles through the slit constraints determined with 
respect to the optical resolution imposed by the analysing apparatus. 

30 The invention may be used in situations requiring non-invasive analysis 

techniques, and may be used in the medical and non-medical fields, in relation 
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to organic and non-organic materials, inanimate or living objects, and crystalline 

and amorphous substances. 

In the medical field, the invention may be used e.g. in the cardiovascular, 

oncological and urological fields in measuring occlusions. 
5 In non-medical fields, the present invention may be used in e.g. the 

characterisation of light metal alloys, e.g. as used in the aerospace industry. 

The invention may be used to measure e.g. stresses and strains in crystalline 

material (the invention analysing the movement of groups of atoms, as opposed 

to detecting individual atoms themselves). 
1 0 The invention may be extended to provide further information about a 

sample e.g. in the y-axis and z-axis. For example, the apparatus could be 

arranged to further scan in the y-direction, and the object or apparatus could be 

rotated to provide other beam paths through the object, and to correlate the 

information obtained in a manner as in tomography. 
15 It is to be understood that various modifications and/or alterations may 

be made to the above without departing from the spirit of the invention as 

outlined herein. 

For example, although discussed in terms of x-rays generated by an x- 
ray tube, any other suitable x-ray source could be used, including e.g. 
20 synchrotron radiation. Also, the monochromator and analyser could take any 
other suitable form, and could e.g. be comprised of double crystals and could 
be made of germanium. Alternative diffraction analysers could also be used 
such as X-ray mirrors. The detector could also take any other suitable form. 



